Forces between chemically structured substrates mediated by critical fluids 
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We consider binary liquid mixtures close to their critical points confined by two parallel, geomet- 
rically flat but chemically structured, substrates. Universal order parameters profiles are calculated 
within mean field theory for periodic patterns of stripes with alternating preferences for the two 
species of the mixture and with different relative positions of the two substrates. From the order 
parameters profiles the effective forces between the two plates are derived. The tuning of Casimir 
amplitudes is discussed. 

PACS numbers: 64.60.Fr, 68.35.Rh, 61.20.-p, 68.35.Bs 



I. INTRODUCTION 

Since chemical structuring of surfaces in the /zm-range 
and below have become feasible, e.g., via self-assembled 
monolayers (SAM) of block-copolymers or mixtures of 
polymers^* 2 ^ micro-contact printing^*^ or exposure 
through a mask of photosensitive surfaces or surfaces 
which change their properties after irradiationj£i£*& 
chemically patterned substrates have gained both prac- 
tical and theoretical interest. Chemically structured 
substrates have applications in micro-reactors, the 
laboratory on the chip, and in chemical sensors^S*ii For 
these systems the local properties of fluids as confined 
by structured walls play an important role. In turn, 
we ^j n gi2 4 i3 i i4 i i5 ag we u ag cr iti ca j adsorption of fluids- 

at chemically structured substrates are in the focus of 
theoretical studies, the latter being the basis for the 
present analysis. While the above studies deal with a 
bulk phase in contact with a single substrate one is also 
interested in thin films confined between two substrates. 
Examples are liquids confined in (slit) pores or thin 
albeit still large on the microscopic scale - wetting films 
for which the vapor phase plays the role of the second 
substrate. For instance there are studies of confined liq- 
uid crystals,' 17 confined block copolymers ^i^S&SL and 
spherical particles with attractive interactions confined 
between two structured substratesi 22 i 23 ' 24 i 25 i 26 i 27 

If a fluid confined between two plates is brought close 
to its critical point, an effective force arises acting on 
the walls due to the boundary conditions imposing a re- 
striction on the spectrum of the critical order parameter 
fluctuations. Since this force has a similar origin as the 
Casimir force between two conducting plates discovered 
by Casimir in 1948, — this force is called critical Casimir 
force and was first predicted by Fisher and de Gennes 29 
and is described by universal scaling functions^SiSiiS For 
fluids confined between two homogeneous substrates pro- 
viding symmetry breaking boundary conditions the uni- 
versal scaling function of the critical Casimir force can be 
calculated analytically within mean field theory^ The 
Casimir forces due to director fluctuations in liquid crys- 



tals confined between a patterned and a homogeneous 
substrate have been recently studied within Gaussian 
approximation^ Here we focus on the universal features 
of the force generated by binary liquid mixtures close to 
their critical point and confined between two chemically 
inhomogeneous substrates with laterally varying prefer- 
ences for the two species of the binary mixture. Within 
mean field theory we calculate numerically the order pa- 
rameter profiles and fully take into account the effect of 
critical adsorption>i§iSi 

The article is organized as follows: In Sec. [n] we recall 
some basic aspects of critical forces between homoge- 
neous substrates in order to set the stage for the inves- 
tigation of the effective forces emerging between inho- 
mogeneous substrates in Sec. 11111 In Subsecs. 1111 Bl and 
IIII CI we study two substrates exhibiting the same and 
opposite patterns, respectively, while in Subsec. lillDl wc 
investigate the forces between one structured and one ho- 
mogeneous substrate. After the comparison of these dif- 
ferent cases in Subsec lill El we conclude in Subsec flll Fl 
with the study of the normal and lateral forces appearing 
between two substrates as function of the misalignment 
of their surface structures. In Sec. llVl we summarize our 
results. 



II. EFFECTIVE FORCES BETWEEN 
HOMOGENEOUS SUBSTRATES 

A. Surface universality classes 

Binary liquid mixtures near their critical point of dcmix- 
ing belong to the so-called Ising universality class which 
also contains one-component (and multi-component) liq- 
uids near their liquid-vapor critical point and uniaxial 
ferromagnets near the Curie point. In this sense in the 
following we refer to binary liquid mixtures without loss 
of generality. Surfaces endow confined systems of a given 
bulk universality class with a fine structure of surface uni- 
versality classes as far as their surface critical behavior is 
concerned. They are characterized within renormaliza- 
tion group theory by the surface enhancement c and the 
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surface field hi. Accordingly, the bulk Ising universality 
class splits up into three surface universality classesii&2ii 
The parameter set (hi — 0, c > 0) characterizes the so- 
called ordinary surface universality class for which the 
order parameter at the surface is suppressed below its 
bulk value. For magnetic systems this is the generic 
case, (hi = 0, c = 0) describes (within mean field the- 
ory) the so-called special surface universality class with 
a multicritical point at which the confined system un- 
dergoes simultaneously a phase transition in the bulk 
and at the surface. For (hi — 0, c < 0) the system be- 
longs to the so-called extraordinary surface universality 
class for which the order parameter is larger at the sur- 
face than in the bulk even above the critical tempera- 
ture for which the bulk order parameter vanishes. This 
is rather uncommon for magnetic systems mirroring the 
choice of the name. In contrast, for fluid systems an 
enhanced order parameter relative to the bulk value is 
the generic case such that this surface universality class 
is then called normal surface universality class which is 
defined by (\hi\ > 0,c = 0). The normal and the ex- 
traordinary surface universality class are equivalent and 
identical at their fixed points, — > oo, c = 0) for fluid 
systems and (hi = 0, c — > — oo) for magnetic systems, 
and thus exhibit the same leading asymptotic critical 
behavior In the following we shall focus on binary 
liquid mixtures close to their bulk critical demixing tem- 
perature T c and exposed to substrates with strong surface 
fields as representatives of the normal surface universal- 
ity class. Our results equally hold for systems at the fixed 
point of the extraordinary surface universality class. 



B. Scaling of the film free energy and the solvation 

force 



The free energy iltot of a binary liquid mixture confined 
between two parallel planar homogeneous substrates of 
macroscopic area A which are separated by a distance 
L and characterized by a surface field hi and h 2 , respec- 
tively, and a surface enhancement c\ and c 2 , respectively, 
consists of a nonsingular, i.e., analytic background con- 
tribution Q ns and a part f2 with a singular dependence 
on temperature: 

fl to t(t,L,hi,h 2 ,ci,c 2 ) 

= Q, na (t,L,hi,h 2 ,ci,c 2 ) + Q(t,L, hi,h 2 ,a,c 2 ) . (1) 

We consider the case t = (T - T c )/T c -> 0* and L > 
atomic length scales for which f2 is expected to exhibit 
universal behavior. Since we focus on the fixed point 
behavior (hi — > ±00, c = 0), the dependences on ci, c 2 , 
hi and h 2 drop out from fl. 

Both contributions to the free energy decompose into four 
distinct contributions: 



U(L) 



LuJb + LJ S> 1 + uj s ,2 + w(L) . 



(2) 



The term ujf, is the bulk free energy density and u s> i and 
u> Sl2 are the surface free energies per surface area A of 
semi- infinite systems bounded by the surface Si and S 2 , 
respectively. The finite-size contribution ui(L) vanishes 
in the limit L — > 00 unless opposing boundary conditions 
create an interface in the middle of the fil m. Conc erning 
this decomposition near T c see, e.g., Refs. 30.31.4ll uib is 
independent of boundary conditions, w s j depends on the 
boundary condition at Si only, and u depends on both 
boundary conditions. Near T c the singular contributions 
exhibit the following scaling behaviors: 



±\d 



a(l - a) (2 - a) 



t 



and 



±\d-l 



a s (l - a s )(2 - a s ) 



(3) 



(4) 



i=l,2, 



where a is the critical exponent of the bulk specific heat 
and a s = a + v is the critical exponent of the surface 
specific heat, £ are the non-universal amplitudes of the 
bulk correlation length t^(t —> ± ) — £o\t\~ v , are 
universal bulk amplitudes, and aj 1 are universal surface 
amplitudes. 

Using the hyper-scaling relation 2 — a — dv^^ the 
singular part of the film free energy can be rewritten as 



1 



k B T c A L d - X I a(l - a) (2 - a) 



a s (l - a s )(2 - a s ) 



(i ± ) d - 1 + e ± (Z ± ) 



(5) 



where & ± (L ) denote the universal scaling functions of 
the finite-size contribution. Since the critical point of the 
film is shifted to a temperature below T c of the bulk, Y 
is an analytic function at T c and thus exhibits the Taylor 
expansion 



Y L 







i=0 



Af(hi = ±00, h 2 = ±oo)(L d 



(6) 



In the case of boundary conditions hi = +00, h 2 = —00 
the analogue of the bulk critical point in the film is 
completely eliminated (see, e.g., Ref. an d references 
therein). Apart from the spatial dimensions and the 
dimensionality of the order parameter the universal 
coefficients depend on the boundary conditions 

applied at both substrates. In Eq. © the coefficients of 
zeroth order 9 ± (0) = A* (hi = ±00, h 2 = ±00) = A++ 
and 9 ± (0) = Aj (hi = ±00, h 2 = Too) = A+~ are 
called Casimir amplitudes so that at T c the finite size 
contribution to the free energy scales as 



<4£ 



L 



-(d-i) 



A. 



T = T r 



A++,A H 



(7) 
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The solvation force F to t(t,L) between the substrates is and 
given by F tot = — ■ Like the free energy it splits up 
into a nonsingular background F ns and a singular part 
F(t, L). The singular part exhibits the scaling behavior 



F(t,L) 



f Xt0 (t,L) = (d-l)L- d f± (L±), (8) 



k B T c A 

with the universal scaling function 
(d- 1)^(1^ 



P ± {w ± -> L ± ,Z ± ) - (L* -w 



± „„±\- 



(15) 



The half-space fixed point scaling function decays expo- 

6 



nentially towards its bulk value P { 
P±( W ± ^oo)-P fc 



where P h + = and P b = 1. 



(16) 



a(l — a) (2 — a) 



D. Mean field theory 



(d-ije^) 



(9) 



In the view of the analysis following later the index 
indicates that homogeneous substrates with no lateral 
structure are considered. One general property of the 
scaling function of the force is that for T ^ T c it 
vanishes exponentially for L — ► oQ-gSLSLMiSl 



f±,o(L — * oo) ~ e . (10) 

At T c the force is determined by the Casimir amplitude: 

/ ± ,o(t = 0,L) = (d-l)iT d A o , (11) 

A = A++A+-. 

The index indicates the case of homogeneous sub- 
strates. 

C. Scaling of the order parameter 

For laterally homogeneous surfaces S± and S2 the or- 
der parameter <p depends on the normal distance z from 
the surface z = and exhibits the following scaling 
properties 



35.36.44.45 



<t>(t,z,L) = a\tfP ± (w ± ,L ± ) 1 w ± = * L± = A 



P : 



(12) 



with P = " P ± and a as the non-universal ampli- 

tude of the bulk order parameter. For increasing separa- 
tion of the substrates the half-space fixed point scaling 
function P^(w ± ) is attained: 



P ± (w ± ,L ± ^^)^P^(w ± ) 



(13) 



For a fixed scaled separation of the two substrates 
the scaling function behaves as: 



P ± {w ± 



0,L d 



(14) 



In the spirit of a systematic field theoretical renormal- 
ization group approach the scaling functions introduced 
above can be determined in lowest order perturbation 
theory, i.e., within mean field theory corresponding to 
e = 4 — d = 0. The scaling variables entering into the 
scaling functions are used with their full scaling form for 
e = l. This approach can be extended also to the case 
of laterally inhomogeneous confining substrates which 
leads to analogous scaling functions with an enlarged 
set of scaling variables. To this end we consider the 
fixed point Hamiltonian TL[<f\ which is a functional of 
the order parameter profile 4>{z). For binary liquid mix- 
tures undergoing demixing it describes the local devia- 
tion of the concentration from the critical concentration 
in the bulk. The dimcnsionless fixed point Hamiltonian 
H[4>] = Hb[4>] + H s [<j>] (providing the statistical weight 
e -WM f or orc [ er parameter <p) separates into the bulk 
part Hb [(f)] confined to the volume V = AL and the sur- 
face part H s [4>] a t the two surfaces with area A&^2& 



4! 



mm = j v & d -\j\z(±{V4>y + T -c 

(17) 

Hs[4>] = f d. d -\ (-MO = 0) - h 2 (f>(z = L)) . (18) 

J A 

Here ry is a vector parallel to the substrates, z indicates 
the direction perpendicular to the substrate, and d de- 
notes the spatial dimension of the system, r is propor- 
tional to the reduced temperature t — T j T ° , the cou- 
pling constant u > 0, which is dimcnsionless for d = 4, 
stabilizes the Hamiltonian %[</>] for temperatures below 
the critical point (T < T c ), and (V</>) 2 penalizes spa- 
tial variations; h± and hi denote the surface fields for 
the substrate at z = and at z = L. Accordingly, the 
Lagrange density Cb[<fi] corresponding to the bulk part 
TLb[<j>] — J dV£b[4>] of the Hamiltonian reads: 



A[0]-i(V0) 2 + I, 



it 
4! 



(19) 



The functional form of the Lagrange density C s [(f)) of the 
surface part Tl s [(f)] is independent of L and leads to the 
boundary conditions, c.f., in Eqs. (|22|l and <|23[) . 
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Within mean field approximation fluctuations of the or- 
der parameter profile are neglected and only the config- 
uration to of the order parameter profile with the largest 
statistical weight e~~ n ^ is taken into account. This mean 
field solution to of the order parameter profile is deter- 
mined by 



sh\ 



= o. 



(20) 



For a system confined between two parallel homogeneous 
substrates Eqs. I|17 )) -Q2u |) lead to a differential equation 
for the mean field profile m{z) of the order parametes 3 ^^ 



-d 2 z m + Tm + — to 3 = 0, (21) 



with boundary conditions 



d z m 



and 



d M z=L = + h 2 ■ 



(22) 



(23) 



E. Calculation of the force from the free energy 

The effective force acting on the two parallel plates is 
the negative derivative of the free energy of the system 
with respect to the distance L between the substrates. 
If this force between the substrates is calculated numeri- 
cally from the difference in the mean field free energies of 
the systems with separation L and L + AL between the 
substrates, numerical problems may arise due to the di- 
verging contributions to the free energy from the surfaces 
for surface fields hip — > ±oo (see, c.f., Eq. (|24ll ) . 
For systems confined between two parallel homogeneous 
substrates with boundary conditions (|22l) and i|23|) the 
free energy of the system per area A and per ^T,. can 
be expressed in terms of the order parameter profiled 



fi[m] 
k B T c A 



h im {0) 

g ( d * m <:„ 



-h 2 m(L) 

L) 2 -^o) 2 



(24) 



4! 



m(z ) 4 



where < zq < L denotes an arbitrarily chosen point 
between the two substrates. Whereas the first two terms 
in Eq. (|24() describe the surface free energy density per 
ks T c , the finite size contribution is contained in the third 
term in Eq. I|24|) . The alternative method for calculating 
the force between the substrates by using the stress tensor 
Tpv (see, c.f., Subsec. IiTF|| yields that the term in square 
brackets in Eq. I|24|l corresponds to the force between the 
substrates^ 3 - In the following we pursue the stress tensor 
approach. 



F. Force calculation using the stress tensor 

This method not only has the advantage that it avoids 
calculating the difference between two diverging free en- 
ergies but also that for determining the force at a distance 
L only the order parameter profile for this distance L is 
needed. Here we present this method for homogeneous 
substrates. It is applicable in the same manner for inho- 
mogeneous substrates (see, c.f., Subsec. Ull A|l . 
The force is a response to the change of the free energy 
of a system if one of the confining substrates is shifted. 
A shift of one of the substrates by an infinitesimal vector 



a(r) 



, /i = 1, . . . , d, leads to a displacement <5r M of 



the components of the vector r, a shift 6m of the order 
parameter, and a change SH in the free energy which is 
given by an integral over the stress tensor and the 
improvement term 



(25) 



where the stress tensor 7^„ and the improvement term 
X^y are given in terms of the Lagrange density C, the 
order parameter to, and their derivatives: 



T 

1 flV 



dC 



j — 



d{d v m) 11 
ld-2 



- SnyC , 



4d- 



(26) 
(27) 



where 8^ u is the Kronecker symbol. 
Now we consider a system confined between two homoge- 
neous substrates and described by the Lagrange density 
in Eq. (|19l) . For this system within mean field theory, 
i.e., (2 = 4, only the component T zz of the stress tensor is 
nonzero: 

T zz = 1 -{d z mf- T -m*-^m\ (28) 

If one of the substrates is shifted by a in normal direc- 
tion (z-direction), a(r) = a6{z — zo)e z , where zq is an 
arbitrary position between the substrates and 9 is the 
Heaviside step function, the Hamiltonian changes by (in 
the laterally homogeneous case = 0) 

SH = J d d r \{d z a z )T zz ] (29) 

d^rn j &z(a8{z-z )T zz ) 

= AaT zz (z = z ) . 

Within mean field theory (<2 = 4), A is a three- 
dimensional surface area. This change SH results in a 
force F± o in z-direction between the two substrates: 



F±,o 



SH 

da 

AT zz (z = z ) . 



(30) 
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where in view of the analysis following below the index 
indicates that homogeneous substrates with no lateral 
structure are considered. For F± t o < (> 0) the confin- 
ing substrates attract (repel) each other. With Eqs. i|28|) 

and (|30j) the force f±,o — k B T° A P er area an( ^ ^bT c is 
given by 



h,o(t,L) = -{d z m{zv)f 



■m(z ) 2 



4! 



m(z ) 4 . (31) 



Equation 131(1 is equivalent to the term in the square 
brackets in Eq. As expected the force f± t o(L) 

exhibits the scaling property described in Eq. JSJ with 

= faltf where £ + = 1/Vt and £" = l/y/2\r\. 
According to Eq. I|31[l the mean field expressions for 
/j_,o and /j_ ; o contain a common prefactor 6/w which 
diverges ~ 1/e for e — * if one inserts for u its fixed 
point value u* ~ e (compare Eq. @ where a ~ e). For 
the quantities studied below we shall consider ratios of 
such scaling functions for which this prefactor drops out 
so that these ratios have finite and well defined mean 
field values in the limit e — * 0. Therefore the mean 
field results for these ratios are expected to provide 
reasonable estimates for their actual counterparts in 
d=3. 

In the case of homogeneous substrates the force between 
the substrates can be calculated directly based on the 
stress tensor without calculating the order parameter 
profile m(z)m& For the fixed point values hi — > +oo, 
/i2 — * +oo, c\ = C2 = the scaling function /jl of 
the force (the index pair ++ denotes that the two surface 
fields are parallel; the symbol ± for T ^ T c is suppressed) 
is given in terms of the complete elliptic integral of the 

' " ' ., with the modulus 



first kind K(k) = ft 
k € [0,1): 



y/l-k 2 l 



(d-l)u 



{2K(k)) 4 k 2 (k 2 - l)-0(-t) 



(L+) 4 



t > t* = - 



where k = k(L + ) is determined implicitly via 



L+(t > f) = Vl(2^(fc)) 2 (2fc2_i)|. 



(32) 



(33) 



The term 



> , ; is the bulk contribution to the force be- 



low the bulk critical temperature. (If one considers the 
slit system together with the thin walls to be immersed 
into the bulk liquid, this bulk contribution to the effec- 
tive force is absent.) For simplicity and practical advan- 
tages below we keep the notation L + also for tempera- 
tures t < 0, with the meaning L + (t < t*) = L~£q /£q~. 




FIG. 1: Universal scaling functions f|J(L) and f±~o{L) of 
the force within mean field theory (see Eqs. (SJ and 13111 ) 
between two parallel geometrically flat and chemically homo- 
geneous substrates at distance L for strong parallel (++) and 
antiparallel (+ — ) surface fields hi and /12. For comparison 
both scaling functions are normalized by the absolute value 
of ff^(0) at T c given by |Aj + | (see Eq. (gUJ). Note that in 
order to avoid breaks in the slope for the curves here we have 
used the same scaling variable L + — L/£ + both above (L > 0) 
and below (L < 0) T c . • and ■ correspond to L/£(t = t*) 
(Eq. lEU) and L/£(t = t°) (Eq. respectively. 



For the temperature range t < t* one has 



t<t\ 
with 



6 



{d-l)u 



(2K(k)) 4 k 2 - 9(-t) 



+ \4 



L+(t < f) = y/(2K(k)) 2 (k 2 + l). 

In the case of strong antiparallel surface fields hi 
ho — ► — 00 one has^ 



(34) 

(35) 
+00, 



f~to(L + ) 

6 



(d- l)u 



(2X(fc)) (1 - fc ) -8(-t) 



t > t° 



with 



L+(t > t°) = ^2{2K{k)) 2 {k 2 + 1) 
For the temperature range t < t° one has 



(36) 



(37) 



fto(L + ) 



G 



{d-l)u 



{2K{k)f~9{-t) 



(L+) 4 



t < t° 



(38) 
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with 

L+(t < t°) = V|2(2^(fc)) 2 (2fc 2 - 1)| . (39) 

Figure n shows these scaling functions. 
For parallel surface fields the force /]^(sign(i)Z + ) is al- 
ways negative, i.e., attractive, and exhibits a minimum 
at a temperature above T c . For antiparallel surface fields 
the force (sign(i)Z + ) is positive, i.e., repulsive. For 

L + — > oo the forces vanish ~ e~^ + = e~ L ^ for t > 
whereas they vanish ~ e~ L £° = er L = er L ^ for 
t < 0. 

At the critical point, the forces are proportional to the 
Casimir amplitudes, which are given by the scaling func- 
tions at L = (see Eq. (IHllh 

/++(L = 0) = A++, (40) 
/£o(£ = 0) = A+-. (41) 

The forces at the critical point are given for k = 1/V2 
and for d = 4 amount to 

< + = -<*(*)) 4 (d^ = - 15 - 7561 ;; < 42 > 

and 

A o"= 1 <^fe)) 4 (^ = 63 " 0242 ^ <«> 

i.e., within mean field theory at the critical point the ab- 
solute value of the force between homogeneous substrates 
with antiparallel infinite surface fields is four times larger 
than for the case with parallel infinite surface fields: 
A+-=4|A++|. 

III. FORCES BETWEEN INHOMOGENEOUS 
SUBSTRATES 

A. Geometry and general expressions 

In this section we consider systems above and at the 
bulk critical temperature. To simplify the notations 
in the following we omit the index + indicating quan- 
tities above the critical point, e.g., L + is replaced by 
L = L/t+ = L/C 

The substrates investigated here exhibit inhomogcncitics 
in one lateral direction (x direction) and translational in- 
variance in the remaining d — 2 directions yi, . . . , yd-2 
(extension H). The lateral substrate inhomogeneities 
form periodic patterns of "positive" stripes (h\ = +oo) 
of scaled width 5+ = 5+/^ and "negative" stripes 
{hi = — oo) of scaled width S- = with periodic 

boundary conditions in x direction (see Fig.[5J). (In the 
preceding paper— these quantities are denoted as S p and 
S n , respectively.) In the following four different combi- 
nations of the substrate structures are investigated: 



(a) y (b) y 




(c) y (d) y 




FIG. 2: Two parallel, geometrically flat, and chemically inho- 
mogeneous substrates at distance L. We consider four generic 
configurations: (a) H — h, (b) +— , (c) periodically structured 
substrate vs. homogeneous substrate, and (d) periodic struc- 
tures shifted relative to each other. On the white stripes 
h\ = — oo, on the grey stripes hi = +oo. 



1. H — h configuration: two identically patterned sub- 
strates with parallel alignment of the infinite sur- 
face fields on the opposing, unshifted stripes. 

2. H — configuration: two identically patterned sub- 
strates with antiparallel alignment of the infinite 
surface fields on the opposing, unshifted stripes. 

3. ph configuration ( "periodic-homogeneous " configu- 
ration): a patterned substrate opposing a homoge- 
neous substrate. 

4. ps configuration ( "periodic- shifted" configuration): 
two substrates with the same periodic pattern of 
stripes of width S but with varying shift D of the 
stripes. 

Various length scales enter describing the inhomo- 
geneities of the substrates, i.e., the set (L,S + ,S-) in 
the cases of the ++ configuration, the H — configuration, 
and the ph configuration. The scaling function f con (L) 
of the force (with con = (++, H — ,ph)) between the two 
patterned substrates depends on the scaled distance L 
and the scaled stripe widths S + = and 5_ = £-/£ 

or equivalently on the set (L, S+/L, S-/S+). The lat- 
ter set exhibits the advantage that upon approaching the 
critical temperature (t — > 0, £ — * oo) only the variable 
L vanishes. The scaling behavior of the normal force 
between two inhomogeneous substrates is given by (com- 
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pare Eq. (jHJl) 

Ff"(f,£,S+,SL) 
k B T c A 



f? n (t,L,S + ,S_) (44) 

(d-i)L- d /r i (i^ + /L,^-/5 + ). 



Without lateral shifts for these configurations 
(++, H — ,ph) the lateral force Fn on is zero. 
For the ps configuration the parameters entering are the 
separation L between the substrates, the stripe width 
S, and the distance D by which the stripes of the second 
substrate are shifted relatively to the first substrate. The 
scaling function f ps (L) of the force is determined by the 
set (L,S/L,S), with S = D/S. In this case the scaling 
behavior of the force normal and lateral to the substrates, 
respectively, is given by (compare Eq. (|5|)): 



and 



F* a (t,L,S,D) 
k B T c A 



Ff(t,L,S,D) 
k B T c A 



= fi s (t,L,S,D) 

= (d-l)L- d fl s (L,S/L,S) 



(45) 



ff(t,L,S,D) (46) 
(d-l)L- d f~P s (L,S/L,6). 



For the configurations H — h, H — and ph the order 
parameter profile exhibits the scaling behavior (with 
con = ++, H — ,ph): 



<t>(t,x,z,L,S+,S-) 



-a\tfP con f v = -,w = -,L = 



z ~ 



(47) 



where P con (v,w, L, S+/L, S-/S+) is a universal scaling 
function and a the non-universal amplitude of the bulk 
order parameter (compare Eq. In the ps case, the 

scaling variables S+/L and S—/S+ are again replaced by 
S/L and 5. 

Quantitative results for the scaling functions of 
the forces follow from the fixed-point Hamiltonian 
H[(f>] =Hb[4>] +H a [<f>] (compare Eqs. JE) and JTSJ). 
Here the integration over the surface area A = N H B 
is Jd^-^n = AT d d - 2 r,| dx, with 5 = 25 for the 
ps configuration and B = 5+ + SL for the other config- 
urations; N is the number of repeat units in x direction. 
The mean field approximation m(x, z) of the order pa- 
rameter profile fulfills the partial differential equation 



(9: 



2 +d 2 )m + Tm + —rr? = . 
' 3! 



(48) 



The accompanying boundary conditions are given by 
Eqs. H22|) and 123|) with h\ and hi varying along x steplike 
at the chemical steps on the substrate. 
From the order parameter profile the force between the 
substrates can be calculated via the stress tensor method 



(see Subsec. Ill F|) . The inhomogeneities in x direc- 
tion lead to further non-vanishing components of the 
stress tensor and of the improvement term (see 
Eqs. D2J), H3) and QZQ). 

A shift of one substrate by a in the z-direction, 
a(r) = a6{z — zo)e z , where zq is an arbitrary position be- 
tween the substrates, yields the normal force f±(L) per 
k B T c and per area A = N B H: 

ML) = ^J dx (T zz (z Q ) + J zz (z )) , 

(49) 

with the components (recall that d Vi m = 0) 



k B T c NBH 



7:\O x m) + —\O z m) — — m — m 

2 K ' V ' 2 4! 



and 



J Z z = \{dxmf + -rndlrn. 



(50) 



(51) 



The lateral force f\\(L) per k B Tc and per area 
A = N B H follows from considering a shift by (3 in 
x-direction, b(r) = f36(z — z )e x : 



F« 



k B T c NBH 
with 



f\\(L) = ~J dx(T xz (z ) + J xz (z )), 

(52) 



%z = {d x m) (d z m) 



(53) 



and 



3 x 



-(d x m) (d z m) - ^md x d z m. (54) 



Within the present mean field level (d = 4) the scaling 
functions f±(L) and f\\(L) as defined in Eqs. (12) - (|4T|) 
can be expressed in terms of the corresponding scaling 
functions P of the order parameter: 



h(L) = 



6 i 4 
— x 



(d- l)u B 

P(w )d 2 v P\ 



(55) 



W I 

1 



3 - - lu-o 2 

; 4 f-B 



(56) 



Here -B = B/^ + and L = denote the scaled exten- 

sions of the system in the x and the z direction, respec- 
tively. 

At the critical point, i.e., for L = the scaling function 
of the force reduces to a generalized Casimir ampli- 
tude A which depends on the type of the inhomogeneities 
of the substrates: 



fx(L = 0) = A<- 



-,ph,ps. (57) 
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B. Forces for the ++ configuration 

In this subsection we consider two identical parallel flat 
substrates with the relative positions of the patterns be- 
ing in phase (see Fig.|^a)). 

First we want to state some general properties of the force 
/]_ ■ The force is symmetric w.r.t. to interchanging 
the width S + of the positive and the width S- of the 
negative stripes: 

/++(£, qi, <fe) = /++ (i, qi q 2 , 1) , (58) 
S+ S- 

This symmetry equally holds in the case of the H — con- 
figuration for the scaling function of the force (see, 
c.f., Subsec. lITlTl . 

If one fixes the width of the positive stripe S+ and con- 
siders the limits S— — » or 5_ — > oo, the force reduces to 
that between two corresponding homogeneous substrates: 

^( x «T'fi:~' ) = I1ML) (59) 

= /++(£) (60) 

This can be seen in Fig. [3J 

At the critical temperature T c , i.e., for L = 0, the scaling 
function of the force depends on S + /L and S-/S+. In 
the subsequent limit of large substrate separations (i.e., 
for S+ /L — ► 0) the dependence on S+ /L drops out so 
that asymptotically the force is governed by a generalized 
Casimir amplitude A ++ (S-/S+) (see, c.f., Fig.EJ). 
Figure shows the dependence of the normal force 
f^ + (L, S+/L, S-/S+) on the scaled distance L between 
the substrates for different ratios S-/S+ and a fixed 
ratio S+/L in comparison with the force f*~^(L) (see 
Fig-Q) for two homogeneous substrates with parallel sur- 
face fields. The strength of the attractive force f^ + (L) 
between patterned substrates is reduced in comparison 
to the strength of the force fxo(L) between homoge- 
neous substrates. At the critical point, i.e., for L = 0, 
the force f^ + (L) between the substrates takes a finite 
negative value, i.e., the force between the substrates is 
attractive. The opposite limit L — > oo can be realized ei- 
ther by increasing the distance L between the substrates 
or by raising the temperature. In both cases the force 
/j^ + (L) between the substrates is expected to vanish ex- 
ponentially (compare Eq. H10|0 which is confirmed by our 
data. For a given distance L the scaled width S+ of the 
positive stripe is fixed (due to the choice S+/L = 1/2 in 
Fig.|3J an d therefore for S-/S+ — ► oo the negative stripes 
become larger and the force between the substrates ap- 
proaches the force /J q (L) = f^(L) between two homo- 
geneous substrates (see Eq. (|tjOH h With increasing stripe 




2 4 6 8 10 

L 



FIG. 3: Scaling function /J + (L) of the force between two 
equally patterned substrates, which are in phase, as a func- 
tion of the scaled distance L = L/£ between the substrates for 
different ratios S-/S+ of the stripe widths (dashed curves) 
and for a fixed ratio S+/L — 1/2 compared with the scal- 
ing function /j_q(-£) of the force between two homogeneous 
substrates (solid curve). The curves are normalized by the 
absolute value of the Casimir amplitude |Aq" + | for the case of 
homogeneous substrates (see Eq. 140t h The surface fields on 
opposing stripes are pointing into the same direction. 



width the order parameter distributions at the chemical 
steps forming the stripes become more independent from 
each other so that the healing of the order parameter pro- 
file becomes more pronounced, i.e., the flattening of the 
cross section at z = const of the order parameter pro- 
file across a chemical step becomes flatter with increas- 
ing distance from the substrate^ (For narrow stripes 
this healing is reduced due to the mutual influence of 
the chemical steps onto each other.) Therefore one can 
conclude that healing of the order parameter distribu- 
tion reduces the absolute value of the attractive force. In 
the limit SL/S+ — » the negative stripes disappear and 
the system again approaches the homogeneous case (see 
Eq. Therefore in Fig.OJthe curve for S-/S+ = 1/3 

lies below that one for S-/S+ = 1/2. For the latter ratio 
the force is weakened most. 

The force f^ + (L) between the substrates also depends 
on the scaled stripe widths 5+ and S_ relative to the 
plate separation L. This is illustrated in Fig. 0] for 
systems with a fixed ratio S-/S+ of the stripe widths 
but with different ratios S+/L. It turns out that for fixed 
distances L > 7 the magnitude of the force increases 
for increasing ratios S+ / L because this implies larger 
stripes. For smaller distances L this trend is again partly 
inverted due to healing effects such that the magnitude of 
the force is no longer smallest for the smallest ratio S+/L. 

The influence of the stripe width ratio S—/S+ becomes 
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FIG. 4: Scaling function /J + (L) of the force between two 
patterned substrates with parallel surface fields normalized by 
|Aj + | for different ratios S+/L and a fixed ratio S-/S+ = 2. 



weaker with increasing scaled distance L between the 
substrates (see Fig. EJ. In the limits S-/S+ — ► 
and S—/S+ — > oo the curves reach the values for 
the force /|q(L) between two homogeneous substrates. 

For L = Fig. [3] renders a universal scaling func- 
tion A, ++ (S—/S+, S+/L) which in the limit of large L, 
S—/S+ — > 0, and 5+ /£ = const or in the limit of large L, 
S—/S+ — * oo, and S+/L — const reduces to the Casimir 
amplitude A^~ + for homogeneous systems. In order to 
obtain the generalized Casimir amplitude A ++ (5_ / 'S+) 
it is necessary to consider the limit S+/L — > which 
is difficult to reach numerically. Approximations thereof 
are shown in Fig. 

This shows how the universal Casimir amplitude can 
be tuned between its homogeneous value Aq" + and its 
maximal value (A ++ ) max ~ via patterning. In the 
limit S+/L — > with S+/S- = const, i.e, for fixed 
stripe widths and increasing film thickness, the periodi- 
cally structured substrates appear as substrates charac- 
terized by laterally averaged, effective homogeneous sur- 
face fields. For S+ = SL the lateral average yields a 
vanishing effective surface field, mimicking a Dirichlet 
(ordinary) boundary condition. For Dirichlet-Dirichlet 
boundary conditions the Casimir force vanishes within 
mean field theory as used here; fluctuations yield a small 
attractive force^S^i This reasoning is consistent with 
the vanishing of A ++ (S-/S+, S+/L) at S+/S- = 1 for 
S+/L -> (Fig. For S+ + S_ the lateral aver- 
age yields a non-vanishing effective surface field of equal 
sign on both substrates, so that, the Casimir amplitude 
remains negative (i.e., attractive) for S+/L — > even 
within mean field theory f Fig. 15)1 . 



FIG. 5: Scaling function (S- / S+) of the force between 
two substrates with their surface field patterns in phase and 
normalized by |Aj + | as a function of the ratio S-/S+ of 
the stripe widths for various fixed distances L between the 
substrates and a ratio S+/L — 1/2. The influence of the ratio 
S-/S+ on the force becomes weaker with increasing distance. 
For S-/S+ = 0, oo one recovers the scaling function for the 
force between two homogeneous substrates, i.e., for each curve 
the values at S-/S+ = and S-/S+ = oo are equal. The 
solid line represents the system at T c , i.e., for L = 0; in this 
case for large L with S+/L = const and S—/S+ — > or 
S-/S+ — > oo one recovers the Casimir amplitude for 
the homogeneous case so that due to the normalization the 
solid curve reaches —1. 




0.5 1 1.5 2 2.5 3 



S-/S+ 

FIG. 6: Generalized Casimir amplitude A ++ (S-/S+, S+/L) 
for a critical fluid confined between two substrates with their 
surface field patterns in phase and normalized by |Aq~ + | as a 
function of the stripe width ratio S-/S+ for different ratios 
S+/L. The curve for S+/L = 1/8 serves as an approximation 
for the actual generalized Casimir amplitude A ++ (£_ / S+) 
corresponding to the limit S+/L — > 0. A ++ (S- / S+) varies 
between Aq" + < for S-/S+ = 0, oo and its maximum value 
(A ++ ) m ax — (see main text). Here and in the following 
figures the curves are smoothly interpolating between the in- 
dicated data points. 
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FIG. 7: Scaling function f± (L) of the force between two op- 
positely patterned substrates (Fig.[5Jb)) normalized by Aq"~ 
as a function of the scaled distance L = L/£ between the sub- 
strates for different ratios S-/S+ of the stripe widths (dashed 
curves) in comparison with the scaling function f*~^{L) of the 
force between two homogeneous substrates (solid curve) with 
a fixed ratio S+/L = 1/2. 



C. Forces for the +— configuration 

In this subsection we study two parallel substrates with 
their chemical patterns being the opposite of each other, 
i.e., a positive stripe of width S+ opposes a negative 
stripe of the same width (see Fig.|2Jb)). 

Figure [7| shows how the scaling function f±~(L) of the 
force in units of A,} - depends on the ratio S-/S+ of the 
stripe widths for a fixed ratio S+/L. As expected this 
force is repulsive. As in the case of the -I — h configura- 
tion the force f±~(L) is reduced compared to the force 
/|q(L) between two homogeneous substrates, takes a 
finite value in the limit L — > 0, and vanishes exponen- 
tially for L — + oo. For S-/S+ — ► oo the force tends to 
the force q(L) between two homogeneous, antagonis- 
tic substrates. For small ratios S-/S+ — > the negative 
stripes vanish as the ratio S+/L is fixed and therefore 
also in this limit the homogeneous case is reached. 
When comparing the forces for systems with the same 
ratio S- / S+ of the stripe widths but with different ra- 
tios S+/L (see Pig- EJ (i- e -> systems with different values 
of scaled widths S + and S- of the stripes) the same be- 
havior as for the ++ configuration is exhibited: for fixed 
distances L > 7 and increasing stripe widths the force 
acting in the system tends monotonicly towards the ho- 
mogeneous regime, for smaller distances L the trend is 
partly inverted due to competing effects caused by re- 
duced healing for very small stripes. 

Figure shows the dependence on the stripe width ratio 
S—/S+. It becomes weaker with increasing scaled dis- 



FIG. 8: Scaling function f^~(L) of the force between two op- 
positely patterned substrates (Fig.|5fb)) normalized by Aq~~ 
for different ratios S+/L and a fixed ratio S-/S+ = 2. 




0.5 1 1.5 2 2.5 3 

S-/S+ 

FIG. 9: Scaling function f±~(S-/S+) of the force between 
two oppositely patterned substrates normalized by Aq as 
a function of the stripe width ratio S-/S+ for various fixed 
scaled distances L = L/£ between the substrates and for a 
ratio S+/L — 1/2. The dependence of the force on the ratio 
S-/S+ becomes weaker with increasing distance. 



tance L between the substrates. The limits S-/S+ — > 
and S-/S+ — > oo correspond to the force f±~^(L) be- 
tween two homogeneous substrates. For L = Fig. |5| 
leads to a universal scaling function A ++ (5_/S+, S+/L) 
which in the limit of large L, S—/S+ — > 0, and 
S+/L = const or in the limit of large L, 5-/5*+ — > oo, 
and S+ / L = const reduces to the Casimir amplitude 
Aq for homogeneous antagonistic systems. A gener- 
alized Casimir amplitude A -1 (S-/S+) as the analogue 
of A ++ (S-/S+) (see Fig. EJ) is obtained in the limit 
S+/L — > 0. A numerical approximation thereof is shown 
in Fig. Like for the ++ configuration (Fig. EJ the 
actual generalized Casimir amplitude A -1 (S-/S+) for 
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FIG. 10: Generalized Casimir amplitude A + ~(S-/S+, S+/L) 
for a critical fluid confined between two oppositely patterned 
substrate normalized by Aq as a function of the stripe 
width ratio S-/S+ for various ratios S+/L. The curve for 
S+ / L — 1/8 serves as an approximation for the actual gener- 
alized Casimir amplitude A + ~ (S- / S+) corresponding to the 
limit S+/L -> 0. A + ~(S-/S+) varies between Aj~ > for 
S-/S+ =0,oo and its minimum value (A +_ ) m i„ ~ (see 
main text). 



the H — configuration vanishes for a stripe width ratio 
S+/S- = 1 in the limit S+/L — > corresponding effec- 
tively to Dirichlet boundary conditions evoking a vanish- 
ing force within mean field theory (Fig. I1U[1 . 



D. Forces for the ph configuration 

The configuration of a patterned substrate at z — par- 
allel to a homogeneous substrate with a positive surface 
field hi = +oo positioned at z = L is of special interest 
because it can be realized by a binary liquid mixture with 
its bulk in the vapor phase forming a wetting film near 
demixing on a patterned substrate at z = with film 
thickness L such that at the free liquid-vapor interface 
the vapor plays the role of the homogeneous substrate. 
This kind of experiment is a natural extension of the crit- 
ical Casimir experiments as carried out in Ref. ^3 on a 
homogeneous substrate. 

Figure ^] illustrates the dependence of the scaling func- 
tion fj^(L) of the force between a patterned and a homo- 
geneous substrate on the stripe widths ratio S-/S+ for 
a fixed ratio S+/ L. For increasing negative stripe widths 
and fixed separations of the plates (i.e., S-/S+ — > oo) 
the force becomes more repulsive and reaches the limit- 
ing case of a system with two homogeneous substrates 
with antiparallel surface fields. In contrast, for vanishing 
negative stripes (i.e., S-/S+ — > 0) the force is attrac- 
tive and converges to the force between two homogeneous 
substrates with parallel surface fields. 



FIG. 11: Scaling function f^_{L) of the force between a pat- 
terned and a homogeneous substrate with a positive surface 
field normalized by |Aj + | as a function of the scaled distance 
between the substrates for different stripe width ratios S-/S+ 
and a fixed ratio S+/L = 1/2. 
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FIG. 12: Generalized Casimir amplitude A ph (S-/S+, S+/L) 
for a critical fluid confined between a patterned and a homo- 
geneous substrate normalized by |Aj + | as a function of the 
stripe width ratio S- / S+ for various ratios S+ / L. As in Fig.|S] 
and llOl the curve for S+/L = 1/8 serves as an approximation 
for the actual generalized Casimir amplitude A ph (S-/S+) 
corresponding to the limit S+/L — > 0. A ph (S~ /S+) varies 



between Aq~ + 



< for S-/S+ = and A J - > for 



S-/S+ — oo. Consequently A ph has a zero and we find 
A ph {S-/S+) = for S- ~ S+. 



At the critical point, i.e., for L = 0, and for large L 
with iS_|_ / L = const the scaling function reduces 
to a universal function A ph (S-/S+, S+/L) (see Fig.[T2J) 
which interpolates between A^~ + at S-/S+ = and 



A 



at 



S-/S+ = oo and therefore has a zero. For 
Sl/L -> 0, according to Fig. d AP h (S-/S + ) vanishes 
for S-/S+ ~ 1. For this stripe width ratio, at the criti- 
cal point the Casimir force decays faster than L~ d . 



12 





-0.2 
-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 
-1.6 



1 

S-/S+ = 1/2 


1 1 1 * 

- t ** « * 
<•**•* ^* 
, *..'* X* 

***••**. X ' 

* * •**" *** S 

y •*' y / ~ 
- - ' ' • ''* • / * / 




"s. 


••*' /* / 

y / s + /l= 

/ 1/8 — 


~ /V 

- ?++ \» 

J±,o 

1 


_y 1/2.-.. _ 

0, oo 

1 I I 



4 6 
L 



10 



The dependence of the force on the scaled distance be- 
tween the substrates is shown in Figs. ED GUI an( i El 
for three different ratios S+/L and three different ratios 
S-/S+. For simultaneously decreasing positive and neg- 
ative stripe widths 5+^0 and 5_ — > the force tends 
- depending on the ratio S-/S+ - towards those for the 
different homogeneous cases, i.e., for S-/S+ > 1 towards 
the force between plates with antiparallel surface fields 
and for S~/S+ < 1 towards the force between plates 
with parallel surface fields. The larger the ratio S+/L, 
the larger is the scaled stripe width S + for L fixed and 
the more developed is the healing of the order parameter 
distribution across the chemical steps. Therefore the ab- 
solute value of the force is reduced as compared to smaller 
stripe widths. 



FIG. 13: Scaling function /^ h (L) of the force between a pat- 
terned and a homogeneous substrate with hi — +oo normal- 
ized by |Aj + | for different ratios S+/L with S-/S+ = 1/2 
fixed. 
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FIG. 14: Scaling function /^ h (Z) of the force between a pat- 
terned and a homogeneous substrate with hi = +oo normal- 
ized by |Ag" + | for different ratios S+/L with S- = S+ fixed. 
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FIG. 15: Scaling function /^ h (Z) of the force between a pat- 
terned and a homogeneous substrate with hi = +oo nor- 
malized by |Aj + | for different ratios S+/L with S-/S+ = 2 
fixed. 



E. Comparison of the configurations 

In this subsection we compare the forces for the afore- 
mentioned configurations H — h, H — and ph. First, we find 
that the forces for the ++ configuration and the H — con- 
figuration are of the same order of magnitude. For the 
same stripe geometries the absolute value of the force 
scaling function f^~(L) for the H — configuration is al- 
ways larger than the absolute value of the force scaling 
function f^[ + {L) for the ++ configuration (see Fig. Htj)>. 
The ratio A+~ (£_/£+, S + /L)/\A++(S-/S+, S+/L)\ of 
the generalized Casimir amplitudes for the ++ configura- 
tion and the H — configuration is smaller than the factor 
4 found for the forces between homogeneous substrates 
(for details see Fig. H7)l . This is caused by the healing 
effect for the order parameter distribution. 
For a ratio S+/L = 1/2 (see Fig. 116( a)) and a ratio 
S-/S+ < 1 the force fj^(L) between the substrates for 
the ph configuration is larger than the force f^ + (L) for 
the ++ configuration but smaller than the force f±~(L) 
for the H — configuration. If the stripes become smaller 
(i.e., S+/L — > 0) (see Fig. IToT bl') the force between the 
substrates for the ph configuration becomes larger than 
for the H — and the H — configuration. 
For the ratio S-/S+ = 1 these effects are even more pro- 
nounced due to the antisymmetry of the order parameter 
profiles: Fig. llftf cl shows the scaling function in the case 
S+/L = 1/2 for which the force between the substrates 
for the ph configuration is even smaller than both the 
force for the ++ and the force for the H — configuration. 
With decreasing stripe widths (see Fig. I16f d)) the force 
between the substrates for the ph configuration is signif- 
icantly larger than for the ++ and the H — configura- 
tion. For these very small stripes and for S-/S+ = 1 the 
structured substrate effectively appears as a surface with 
a surface field tending to zero which is effectively equiva- 
lent to the Dirichlct boundary condition. In contrast the 
homogeneous substrate in the ph configuration still has 
a positive surface field such that for the ph configuration 
a force remains even in the limit S+ / L — > 0. 
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L L 

FIG. 16: Comparison of the absolute values of the scaling functions f± n (L) for the forces between patterned substrates (where 
con = ++,H — ,ph) normalized by |A ( |+|, (a) for the ratios S+/L = 1/2 and S-/S+ = 2, (b) S+/L = 1/8 and S-/S+ = 2, (c) 
S+/L — 1/2 and S- = S+, and (d) S+/L = 1/8 and S- = S+; on this scale the values for both the ++ and the H — configuration 
cannot be distinguished from zero. 



The Casimir amplitude A° + of the force between a sub- 
strate with ordinary (Dirichlet) boundary condition and 
a substrate with + boundary condition within mean field 
theory is given by2£ 

A°+ = -jA++. (61) 

The value /f (0)/|A++| for S+/L = 1/8 in Fig.QHd) is 
still above this value because the limit S+/L — > is not 
yet reached. 

Figure El shows the dependence of the ratio 
A + -(S-/S + ,S+/L)/\A ++ (S-/S + ,S+/L)\ of the 
generalized Casimir amplitudes for the H — con- 
figuration and the H — \- configuration on the stripe 
width ratio S-/S + . For S-/S+ = 0, oo the sub- 
strates are homogeneous and therefore in these 
limits A+-(S-/S + ,S + /L)/\A++{S-/S + ,S+/L)\ = 
A^/IA^+I = 4. (Since both A++ and A+" are very 



small for S-/S+ ~ 1 (see Figs. IdI and 1100 . the numerical 
estimates for this ratio become unreliable there.) 



F. Lateral shift 

Finally we consider two parallel substrates endowed 
with the same pattern but laterally shifted relative to 
each other by a distance D in x direction, D = D/£ 
(see Fig. Hfd)). In order to limit the number of rele- 
vant length scales here we choose S + = 5_ = S and 
S+ = S- = S, respectively, so that with the relative 
shift 6 = D/S = D/S the set of variables is (L, S, D) or 
(L,S,5). This set of variables is chosen to simplify the 
presentation of the scaling functions of the force but it 
has the disadvantage that upon approaching the critical 
point both the scaled distance L and the scaled stripe 
width S vanish which makes it difficult to calculate the 
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FIG. 17: Ratio A + -(S-/S' + ,S' + /L)/|A ++ (S'-/5' + ,S'+/L)| of 
the generalized Casimir amplitudes for the +— configura- 
tion and the ++ configuration as a function of the ratio 
S-/S+ of the stripe widths. The limits S-/S+ — 0, 00 corre- 
spond to homogeneous substrates for which the ratio equals 
Aj"/|A^ + | = 4 (see Subsec. UTTl . It is difficult to evalu- 
ate numerically this ratio for S-/S+ — 1 because there both 
A+-(S-/S +y S+/L) and A++(S-/S+, S+/L) are very small. 



forces between the substrates at the critical point. In 
order to obtain the forces at T c the set of variables 
(L, S/L, D/S) = (L, S/L, S) in which only the scaled dis- 
tance L vanishes in the limit T — > T c is chosen. For this 
configuration we study the force f±(L, S,S) normal to 
the substrate like for the previous configurations. But in 
this case in addition a lateral force fn S (L,S,5) parallel 
to the substrates emerges. 



1. Normal force 

Figure ITSl illustrates the influence of the relative shift on 
the normal force f^(L, S, 5). For a vanishing shift 5 = 
the system is in the patterned H — h configuration and the 
force is attractive for all scaled distances L between the 
substrates. With increasing shift S (i.e., D — > S) the force 
increases and becomes repulsive for all scaled distances 
L until the corresponding patterned H — configuration is 
reached for 5=1. For 6 > 0.4 the force becomes repulsive 
for all L. 

The influence of the scaled stripe width S on the scal- 
ing function f^{L, S, 5) is shown in Fig. 1191 through the 
dependence of on the ratio L/S. The curves for the 
scaled stripe widths S = 1.0, S = 2.0, and S = 3.0 with 
the same relative shift 5 = 0.2 do not coincide. This 
demonstrates that the force f±(L, S, 5) does not reduce 
to an effective scaling function of only the two variables 
L/S and 6. 

In Fig. 1201 the data contained in Fig. ^| are shown as to 



FIG. 18: Scaling function (L, S , 5) of the normal force 
between two substrates patterned identically with positive 
and negative stripes of widths S+ — S- — S and S = 2.0 for 
various lateral shifts S = D/S = D/S. The scaling function 
of the force is normalized by |Aq" + |. The values for L = are 
obtained by taking the values for small L and extrapolating 
them linearly to L = 0. 
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L/S 

FIG. 19: Scaling function of the normal force between 
two substrates patterned with positive and negative stripes 
of equal width S — 1.0 (solid line), S = 2.0 (dashed line), and 
S = 3.0 (dotted line) as function of L/S and with a relative 
shift of S = 0.2. Here is considered as a scaling function 
of the scaling variables L/S, S, and S. The scaling function 
of the force is normalized by |Aq" + |. 



visualize the dependence on the relative shift S = D/S. 
For small scaled distances L the force scaling function 
depends strongly on the relative shift, with the curves in- 
terpolating between the value of the force scaling function 
f~^ + (L) for the patterned ++ configuration [8 = 0) and 
of the force scaling function fj~(lS) for the patterned 
H — configuration (5 = 1). The dependence on 6 disap- 
pears with increasing scaled distance L — > 00 between 
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FIG. 20: Same data as in Fig. 1181 shown as function of the 
relative shift 8 = D/S. 



the substrates. 

The normal force at the critical point is characterized by 
a generalized Casimir amplitude A ps (5, S/L) shown in 
Fig. for the ratios S/L = 1/2, 1/4, 1/8. The general- 
ized Casimir amplitude A ps (5,S/L) varies as a function 
of the relative shift S between A++(S-/S+ = 1, S+/L) 
for S = and A+-(S-/S+ = 1,S+/L) for 5 = 1. Like 
the generalized Casimir amplitudes A ++ (S-/S + , S + /L) 
and (S-/S+, S+/ L) vanishing for a stripe width ra- 
tio S-/S+ = 1 (Figs. El and Uni) in the limit S+/L -> 
- resembling effective Dirichlet boundary conditions - , 
the generalized Casimir amplitude A ps (S, S/L) vanishes, 
within mean field theory, in the limit S/L — > 0, too, for 
all values of the shift S (recall that here S+ = SL = S). 



2. Lateral force 

The scaling function f ps (L, S,S) of the lateral forces is 
calculated according to Eqs. (|S"2|l and lj5U|) . Figure 
shows the dependence of the scaling function f^ s (L, S, S) 

for the lateral force on the scaled distance L between the 
substrates for a fixed scaled stripe width S = 2.0. For 
the relative shifts 5 = 0, 2 and 6=1, which correspond 
to the patterned H — h and the patterned H — configura- 
tion, respectively, the lateral force f ps (L,S,S) vanishes 
identically. With increasing shifts up to S « 0.5 the force 
first becomes increasingly restoring and for shifts S > 0.5 
it decreases again towards 0. For shifts < S < 1 the 
direction of the force is opposed to the direction of the 
shift, whereas for shifts 1 < 6 < 2 the force acts into the 
direction of the shift . 

At the critical point, i.e., for L — and S — ► oo in order 
to keep S — 2 as in Fig.[23the patterning becomes irrele- 
vant such that the lateral force attains 0. With increasing 
distance L the force reaches a minimum (maximum) for 
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FIG. 21: Generalized Casimir amplitude A ps (S, S/L) for a 
critical fluid confined between two identically structured sub- 
strates with relative shift 8 normalized by \Ag + \ as a function 
of the relative shift S (a) for a ratio S/L = 1/2, (b) for a ratio 
S/L = 1/4, (c) for a ratio S/L = 1/8; S+ = S- = S. 



0<(5<l(l<<5<2)ata certain L(S) > and in the 
limit L — > oo the force vanishes exponentially due to the 
same arguments which apply to the normal forces. Thus 
as a function of L the lateral force reflects the behavior 
of the normal force shown in Fig. ^3J 
The dependence of the scaling function of the lateral 
force f ps (L,S,5) on the relative shift 8 is illustrated 
directly in Fig. [23 for various fixed scaled distances 
L between the substrates and for S = 2. The force 
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FIG. 22: Scaling function ff(L,S,8) of the lateral force, 
normalized by |Aj + |, between two identically patterned sub- 
strates as function of the distance L between the substrates for 
various relative shifts S and for a scaled stripe width S = 2.0. 



FIG. 23: Scaling function fF" (L, S , 8) of the lateral force 

normalized by |Aj + | between two identically patterned sub- 
strates as a function of the relative shift 8 for different scaled 
distances L between the substrates and for a scaled stripe 
width S = 2.0. 
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attains a minimum for S 0.5, i.e., for this value the 
restoring force is strongest. For 8 ss 1.5 the repulsive 
force attains a maximum. For L <^ S the force scaling 
function exhibit plateaus as function of 8 around <5 = 0.5 
and (5 = 1. 5. At 8 = 1 the force scaling function has 
an inflection point which corresponds to an unstable 
configuration. 

The singular free energy cost for shifting the two sub- 
strates relative to each other is given by (see Eq. 147(0 



Uf(D, L, S, t) = k B T c ^^-(d- l)Uf(8, L, S) , (62) 

where N is the number of repeat units of width 2S and 
H is the extension of the system in the translational in- 
variant direction so that A — 2NSH. The dimensionless 
scaling function 



U* S (8,L,S) = - / d6'f* a (8',L,S), 



(63) 



is shown in Fig. [21 

For given lateral features N, S, and H one has 
2NS 2 H/L d < 1 for L sufficiently large so that Uf < 
ks T c which implies that the structured substrates can 
float against each other via thermal activation if they 
are not externally fixed. For N — 50 repeat units con- 
sisting of a positive and a negative stripe each of width 
S = 100 nm and for an extension H = 10 fmx the above 
condition is met for L > 2 fim. 

The scaling function U^ s (8, L, S) of the free energy ex- 
hibits a maximum U max for <5 = 1, i.e., D = S. This 
maximum is the activation barrier which one has to over- 
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FIG. 24: The scaling function U(S,L, S), normalized by 
|Aq" + |, of the singular contribution to the cost of free en- 
ergy to shear the patterned substrates relative to each other 
by the relative shift S = D/S. The scaling function has its 
maximum at S = 1, minima are at S = and 8 = 2, and 
turning points at 8 = 0.5 and 8 = 1.5. 



come to shear the patterned system from the ++ config- 
uration at 8 — over the H — configuration at 8 = 1 to 
the ++ configuration at (5 — 2. Here the absolute values 
of the free energy do not matter. The turning points at 
8 = 0.5 and 8 = 1.5 of the scaling function of the free 
energy correspond to the maximal lateral force between 
the substrates. 

Figure 1251 shows the universal activation barrier 
U m ax{L, S) with its dependence on the scaled distance 
L between the substrates for a fixed scaled stripe width 
5 = 2. This scaling function vanishes for L — > and 
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FIG. 25: Scaling function Umax{L, S), normalized by |Aj + |, 
of the activation barrier as a function of the scaled distance 
L for a fixed scaled stripe width S = 2.0. 



FIG. 26: Scaling function /jf iw of the lateral van der Waals 
force as a function of the relative shift 8 for L — S (see 
Eqs. (I5H and jgg)). 



L — * oo and, for this choice of 5, exhibits a maximum for 
L«2. 



3. Comparison with lateral van der Waals forces 

As described in the Introduction the chemical pattern- 
ing requires at least one monolayer of the substrate to 
be composed of alternating stripes of distinct chemical 
species which provide the contrast in preference for the 
two types of particles forming the binary liquid mixture 
exposed to this lateral structure. Since these two sub- 
strate species interact not only with the fluid particles 
but also among each other, their inhomogeneous lateral 
distribution gives rise to a lateral force even in the ab- 
sence of the fluid between the two confining substrates. 
This lateral force stems form the direct van der Waals or 
dispersion forces between the substrate particles and pro- 
vides a nonsingular background contribution which adds 
to the singular lateral force discussed above. This non- 
singular contribution can be determined separately by 
taking the fluid out of the system. 

In order to estimate the background force we consider 
two types of particles i,j = +, — interacting pairwise 
via Lennard-Jones potentials which are characterized by 
energy and length parameters ey and cr^ , respectively. 
Their actual number density on the stripes are rji. Since 
L <Jij , only the attractive part of the Lennard-Jones 
potentials matters. As discussed in Appendix ^] (see 
Eq. (|A18() ^. the lateral van der Waals force (for d = 3) 
has the form 



that F^ dw vanishes if the chemical contrast disappears, 
i.e., for + = — ; fn is a dimensionless scaling function 
(see Eq. f£T?ty): 



fvdW ( L A _ 157T \ - 

J " S' 16 ^ 

x / 771— — OO 

dCi / d( 2 

'0 JS+2m 



(65) 



C1-C2 



(§) +(Ci-C 2 ) 5 



For L/S = 1 the scaling function fX dW (see Fig.EHJ at- 
tains maximal values of the order of 0.3. We estimate 
the prefactor E/S 5 as follows. The interaction energy 
is typically in the order of e/fcs = 300if (see Table I in 
Ref. lifll) . the length <Xy is comparable with the atomic 
diameters of the interacting particles, cry = 0.5 nm, and 
in the case of a 2-dimensional closed packing the areal 
number density rji of the particles is rji = v3/(6-R 2 ), 
with the radius Ri — cru/2 of the particles. Assum- 
ing stripe widths of S — 100 nm one has typical values 
E i:j /S 5 ~ 5.5 • 10- 4 pN/(^m) 2 . Thus the lateral van der 
Waals force is of the order of 1.7 • 10~ 4 pN/(^m) 2 . The 
actual value is smaller due to the differences of E^j form- 
ing E. 

According to Eq. I|47() the singular lateral force is given 
by 



Ff(t,L,S,D) _ ( d -l)k B T c 



2NSH 



ps 



L, S T ,S 



(66) 



F v ' 



2NSH 



(D,S,L) 



E - 



vdW 



S' 



(64) 



where E = E. 



+^ 



-2E+-+E — and E itj = 46^^^^, so 



For a ratio S/L = 1 the scaling function /IP 8 attains 

maximal values of the order of | Aj + 1 (see Fig. I23|l with 
|A++| = 0.32^ For d = 3, T c = 300 K, and a separation 
L = 100 nm of the substrates, the critical force is of the 
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order of 1.3 pN/(^tm) 2 . According to this estimate for 
ratios S/L ~ 1 the direct lateral van der Waals force due 
to the chemical patterning is negligible compared with 
the critical structural force. 



IV. SUMMARY 

We have studied the singular contributions to the ef- 
fective forces acting on chemically inhomogeneous sub- 
strates which confine a binary liquid mixture near a con- 
tinuous demixing transition. Geometrically flat and par- 
allel substrates with periodic chemical patterns of pos- 
itive stripes (surface field hi — > +oo) of width S+ and 
negative stripes (hi — > — oo) of width S- separated by a 
distance L are considered. 

Four basic configurations of the two substrates are stud- 
ied. First, two substrates exhibiting the same stripe 
patterns, i.e., a positive stripe opposite of a positive 
stripe, which we call "++ configuration", are considered 
(Fig- Eta) j Subsec. IlIIBJ) . Second, two substrates with 
opposite stripe patterns, i.e., a positive stripe opposed to 
a negative one, called the "H — configuration", are studied 
(Fig.|3b), Subsec. lITIT^ . This is followed by investiga- 
tions of the forces between a structured and a homoge- 
neous substrate, called the u ph configuration" (Fig.EJc), 
Subsec. IIII Dfl . For two substrates with the same but 
misaligned stripe patterns, called the ll ps configuration" 
(Fig- Hid), Subsec. Illl the normal and lateral forces 
depend on the misalignment of the surface structures. 
The universal behavior of the order parameter profiles 
and of the effective forces acting between the substrates 
is determined by universal scaling functions (Eqs. (|44|l - 
1)47(1 ). In order to obtain explicit results for these scaling 
functions the order parameter profiles are calculated 
numerically within mean field theory from which the 
forces between the plates are derived via the stress 
tensor (Eqs. O - © - ©). 

The scaling function of the force between two periodically 
structured substrates in the ++ configuration depends 
on the scaled distance L = L/£ between the substrates, 
where £ is the bulk correlation length, as well as on the 
ratio S-/S+ of the stripe widths (Figs. and [SJ and 
the ratio S+/L (Fig. 0J . In order to obtain well-defined 
mean field values for the spatial dimension d = 4 the scal- 
ing functions of the force between structured substrates 
are expressed in units of the universal Casimir amplitude 
A,| + characterizing the scaling function of the force be- 
tween homogeneous substrates with parallel surface fields 
(Fig. [l]and Eq. igOJl). In the limit S-/S+ -> and 
S—/S+ — * oo with fixed ratio S+/L the scaling function 
of the patterned substrates approaches the scaling func- 
tion of the homogeneous substrates. This limiting case is 
also reached for an increasing ratio S+/L. 
At the critical point the force is determined by a 
universal scaling function of the scaling variables S- / S+ 
and S+/L, which reduces to the Casimir amplitude Aq + 



in the limits S-/S+ -> and S-/S+ -> oo (Fig.©. The 
limit L — * oo of this scaling function can be interpreted 
as a generalized Casimir amplitude A ++ (S-/S+). Thus 
the chemical patterning of the substrates allows one to 
tune the universal Casimir amplitude within a finite 
range of negative values (Fig. © . 

The force between two periodically structured substrates 
in the H — configuration is governed by a scaling func- 
tion which depends on the same scaling variables L, 
S-/S+ (Figs. □ and ©, and S+/L (Fig. EJ. Like for 
the H — h configuration it is reduced compared to the 
scaling function of the force between two homogeneous 
substrates with antiparallcl surface fields. Analogously 
to the ++ configuration, the scaling function of the 
H — configuration attains the scaling function of the 
homogeneous substrates in the limits 5_ / S+ — > and 
S-/S+ — > oo with a fixed ratio S+/L and in the limit 
S+ [L — > oo for a fixed ratio 5_ / S+ . 

As in the case of the ++ configuration, a generalized 
Casimir amplitude A A (S-/S+) can be defined also for 
the H — configuration. This follows from the universal 
scaling function describing the force at T c which depends 
on the scaling variables S—/S+ and S+/L. In the limits 
S-/S+ — > and S-/S+ — > oo it reduces to the Casimir 
amplitude Aq" > (Eq. (^lfl ) characterizing the 
force between homogeneous substrates with antiparallel 
surface fields (Fig. I1U[1 . Here the chemical patterning 
allows to tune the universal Casimir amplitude within 
a finite range of positive values (Fig. ITU)) . The ratio 
of the generalized Casimir amplitudes for the H — and 
the ++ configuration as function of S-/S+ is shown in 
Fig. El 

The ph configuration of a critical fluid confined between 
a periodically patterned and a homogeneous substrate 
with positive surface field interpolates, as a function of 
the ratio S-/S+ of the stripe widths, between the case 
of two homogeneous substrates with parallel surface 
fields (S-/S+ — > 0) and the case of two homogeneous 
substrates with antiparallcl surface fields (S-/S+ — > oo) 
(Fig. El Figs. IT31- IT5|) . Accordingly the corresponding- 
generalized Casimir amplitude A ph (S- / S+) interpolates 
between the Casimir amplitude Ag + < at S-/S+ = 
and the Casimir amplitude Aq~ > at S-/S+ — > oo. 
Figure El shows universal scaling functions at T c for 
non- vanishing ratios S+/L, which converge to the 
generalized Casimir amplitude A ph (S-/S+) in the limit 
S+/L — > 0. In this case the chemical patterning allows 
one to tune even the sign of the Casimir amplitude and 
to drive the leading critical Casimir force to zero. This 
opens the possibility to observe the leading correction 
term which is typically masked by the dominant term. 

The absolute value of the scaling function of the force 
between substrates in the H — h configuration is smaller 
than the corresponding one for the H — configuration 
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(Fig. In the case of very narrow stripes, i.e., for 

S+/L <C 1 fFigs. Uffl bl and (d)) the absolute value of the 
scaling function of the ph configuration is significantly 
larger than the one for the ++ configuration and the 
one for the H — configuration. The reason for this is that 
the surface fields of the structured substrates effectively 
tend to zero for decreasing stripe widths - mimicking 
Dirichlet boundary conditions - whereas the surface field 
of the homogeneous substrate of the ph configuration 
remains finite. 

If the structures of the patterned substrates are not in 
phase but shifted relative to each other, an additional 
lateral force emerges which acts as to restore a configu- 
ration in which the structures are again in phase. In this 
ps configuration the scaling functions of the force depend 
additionally on the relative shift 5 = D/S of the stripes 
of one substrate with respect to the stripes of the other 
substrate (Fig.^d)). In order to limit the number of the 
relevant length scales here we have chosen S- = S+ = S. 
The scaling function of the normal force as a function of 
the relative shift S transfers the scaling function for the 
H — h configuration for 5 = into the scaling function for 
the H — configuration for 5 = 1 fFigs.ir%land l2U|) . 
The importance of the relative shift 5 decreases with in- 
creasing scaled distance L between the substrates. The 
influence of the scaled stripe width S = S/£ is visible 
through the dependence of the scaling function on the 
ratio L/S (Fig. I19|) with the conclusion that the scaling 
function does not reduce to an effective scaling function 
of only the two variables L/S and S. 
The generalized Casimir amplitude A ps (S/L, 6), char- 
acterizing the scaling function of the normal force at 
the critical temperature fFig. I21fl . vanishes in the limit 
S/L — ► 0. For positive and negative stripes of the same 
width, in this limit the surface field effectively tends 
to zero corresponding to Dirichlet (ordinary) boundary 
conditions for which the force vanishes within mean field 
theory. 

Figures H3 El El an d |^ demonstrate the tunability 
of the Casimir amplitudes via the chemical patterning 
of the substrates. The fact that within mean field 
theory the Casimir amplitudes A ++ (S-/S+, S+/L), 
A + -{S-/S+,S+/L), and Ap s (S,S/L) vanish for 
S+/S- = 1 in the limit S+/L -> and S/L -> 0, 
respectively, resembling Dirichlet boundary conditions, 
is remarkable. In general in classical fluid systems 
the substrate potential breaks the symmetry of the 
order parameter such that these systems are unable 
to exhibit Dirichlet boundary conditions. Only in the 
case of quantum fluids, i.e., in critical superfluid He 4 
films, Dirichlet boundary conditions show xvpM^ The 
results shown in Figs.[fjlEIEl and 12 II visualize that for 
suitably chosen chemical patterns Dirichlet boundary 
conditions effectively emerge even for classical fluids. 

The scaling function of the lateral force vanishes for rela- 



tive shifts 6 = 0, S = 1, and 8 = 2 (Figs. El and El- For 
shifts < 8 < 1 the scaling function of the lateral force is 
negative, i.e., restoring. For 1 < 8 < 2 it is positive, i.e., 
the force acts into the direction of the shift. For 8=1 the 
scaling function exhibits an inflection point, which means 
that the corresponding H — configuration is an unstable 
configuration. 

The free energy required to shift the patterned substrates 
relative to each other (Eq. 16211 1 is given by a scaling func- 
tion (Eq. (|63[) and Fig. 12411 which shows a maximum at 
D ~ S and exhibits inflection points at D/S = 0.5 and 
D/S = 1.5 corresponding to the occurrence of maximal 
forces. The maximum of this free energy, i.e., the acti- 
vation energy needed to move a system from a H — h con- 
figuration to the adjacent H — configuration via shifting 
the substrates by D = S, depends on the distance be- 
tween the substrates and vanishes for increasing distances 
(Fig. El- 

The background contribution from direct van der Waals 
forces (Fig. I26|) . generated by the chemical patterning 
within a monolayer covering the substrates, is negligi- 
ble for ratios S/L ~ 1 compared with the critical lateral 
force. 



APPENDIX A: LATERAL VAN DER WAALS 
FORCE 



The lateral van der Waals force F, 



vdW 



dU?, d 



-SB- 18 g^- 



erated by the direct interaction between the two sub- 
strate species forming the chemical pattern. Here we as- 
sume that this chemical pattern is confined to a mono- 
layer and that the corresponding pair potentials are given 
by Lennard- Jones potentials: 



U, 



4e, 



12 



i,je {+,-}. (Al) 



Since L ^S> for the lateral van der Waals force only 
the attractive part of the Lennard-Jones potentials is im- 
portant. Assuming pairwise additivity the leading con- 
tribution to Un dw is given byi^i 



U 



vdW 



dxidyi / dx 2 dj/ 2 
Si J s 2 

E(x\,X2,D) 

[L 2 + {xi - x 2 y + (yi - y 2 ) 2 } 3 ' 



(A2) 



with x\ and x 2 denoting the direction of the inhomo- 
geneity on the substrate Si and S 2 , respectively, and yi 
and y 2 the directions perpendicular to it (see Fig. EJ. 
Depending on the lateral position, E{x\, x 2 , D) equals 
E++, E+-, or (see, c.f., Eqs. {SJ and (X8j|). The 
contributions are given by Eij = 4ejjcr?-%77 J - ! where 
r\i denote the areal number densities of the particles. 



In order to carry out the quadruple integral in Eq. (|A2I) 
we first consider the interaction Uu g of a positive stripe 
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on the substrate Si of width 5 and extension H in the 
2/1 direction with the whole substrate 62 with infinite 
extensions in the x 2 and y 2 directions: 
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dxi / dyi / dx 2 / dy 2 

Jo J-oo J-00 

E(xi,x 2 ,D) 

[L 2 + (xi - x 2 ) 2 + (j/i - y 2 ) 2 f ' 



(A3) 



With 



+00 



dx- 



(a 2 +x 2 ) 3 8a 



S (A4) 



Eq. 1A3H gives 



vdW 



— H I dx\ I dx 



E(x 1 ,x 2 ,D) 
[L 2 + (xi - x 2 ) 2 ]i 



(A5) 



The integration over the x 2 direction can be split up 
into integrations over the stripes of the substrates S 2l 
for which the strength E(x±, x 2l D) of the interaction is 
constant: 

xi G [0 , 5] , x 2 e [D + 2m5 , D + (2m + 1)5] : 
E(xi,x 2 ,D) = E ++ 

x x G [0 , 5] , x 2 € [D + (2m + 1)5 , D + (2m + 2)5] : 
E(x 1 ,x 2 ,D) = E+_, (A6) 

where m e Z. Thus one has for Uf™ ( Eq . 
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l) JD + (2m+l)S [L 2 + (Xi — X 2 ) 2 Y- 



(A7) 



Now we consider the interaction energy Uu g of a nega- 
tive stripe on the substrate 5i of width 5 and extension 
H in the 2/1 direction with the whole substrate 52 with 
infinite extensions in x 2 and 2/2 direction. The strength 
of the interaction is given by 

an G [5 , 25] , x 2 G [D + 2mS , D + (2m + 1)5] : 
E(xi,X2,D) = E+- 

x-l G [5 , 25] , £2 G [D + (2m + 1)5 , D + (2m + 2)5] : 
E(x 1 ,x 2 ,D) = E—, (A8) 



such that J/j!"^ is described by: 
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(A9) 

The sum of the interaction energies Ung (Eq. (IA7jO and 
( E< 1- ®) is 

r D+(2m+l)S 



m=—oo 

dx 2 I / dxi 



D+2mS 



E 



+ + 



2S 



[L 2 + ( Xl - x 2 ) 2 }-2 

E, 
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E, 



dx 2 / dxi — 

D+(2ra+l)S WO [L 2 + (Xx ~ X 2 ) 2 ] 



2S 



dxx 



E. 



s [L 2 + ( Xl - x 2 ) 2 



3?r 



B E 



m— — 00 
D+(2m+l)S 



E ++ - E + _ 



dx 2 / dxi 

D+2mS JO [L 2 + (Xx - X 2 ) 2 ] 5 



E. 



+- 



D+(2m+l)5 <-2S 

dx 2 / dsi — 

D+2mS JO [L 2 + (XX — a; 2 ) 2 ] 2 

n + (2m+2)S |-2S g _ £^ 



dxo, I dx^ 



D+(2m + X)S JS [L 2 + (Xx - X 2 ) 2 } 

D+{2m+2)S r2S j£ 

dx 2 / da?i 

D+(2m+l)S Jo [L 2 + (xx - a:2) 2 ] 
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The transformation X\ = xx — 5+ and A2 = x 2 — 5+ in 
the third term renders: 



,S- 



rrvdW , jrvdW _ ^7T 

u \\,s++ u \ 

r D+(2m+X)S ,- s 



* E 



m=— 00 

i , , E ++ - 2E + - + E- 
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D+2mS Jo [L 2 + (xx - X 2 ) 2 Y- 



E, 



D+(2m+2)S r2S 

dx 2 / dxx — 

D+2mS Jo [L 2 + (XX — X2) 2 ]'- 



(AH) 
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Scaling the lengths with the stripe width S, 



The sum over the second integral reduces to: 



leads to 



S ' ^ 2 5 ' S" 



(A12) 



u \\,s+ + u \\,s. 



37T_ff 

" ft 53 2^ 



00 4+(2m+2) 
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E 4 
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d( 2 / dCi 



£++ - 2E, + E— 
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lim 
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(A14) 
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S+(2m+2) <-2 

dC 2 / dCi 

S+2m JO 



E, 



(I) +(d-e.) s 



(A13) 



J 



The interaction energy UV w /N of a positive and a negative stripe on one substrate, where N is the number of repeat 
units, with the structure on the other substrate per area 2SH reads 

fjvdW 



2NSH 16 S 4 

3tt E. 



Tfi ^ E / d ^ / d( ^ 2 



rn— — 00 



1 



J S+2m 



16 S 4 
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dC 2 



J S-2m 



#J + (Ci - C2) 2 
1 

(#) 2 + (Ci-c 2 ) ; 



(A15) 



where E = E ++ - 2E + _ + 

The lateral van der Waals force is determined by the derivative of UV with respect to the shift D: 
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3tt E 
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rn— — 00 



d_ f 1 
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dCi 
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(A16) 



The last term in Eq. I|A16J| vanishes: 



lim / dCi 

'0 J5-2r. 



<5+(2m+2) 
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# +(Ci-C 2 



- =0. 
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Therefore the lateral van der Waals force Fn per area 2NSH between two structured monolayers is given by 



2NSH S 5 



rdW 



(A18) 
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with the scaling function 
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15tt 
"16" 



E 



m— — og 



1 /.(5+2m+l 

dCi / ckV 

J S+2m 



C1-C2 



(§) +(Ci~C 2 ) 5 



(A19) 



1 F. S. Bates, Science 251, 898 (1991). 

2 G. Krausch and R. Magerle, Adv. Mat. 14, 1579 (2002). 

3 M. Boltau, S. Walheim, J. Mlynek, G. Krausch, and U. 
Steiner, Nature 391, 877 (1998). 

4 A. Kumar and G.M. Whitesides, Appl. Phys. Lett. 63, 
2002 (1993). 

5 Y. Xia, J. A. Rogers, K.E. Paul, and G.M. Whitesides, 
Chem. Rev. 99, 1823 (1999). 

6 B. Zhao, J. S. Moore, and D. J. Beebe, Anal. Chem. 74, 
4259 (2002). 

7 R. Wang, K. Hashimoto, A. Fujishima, M. Chikuni, E. 
Kojima, A. Kitamura, M. Shimohigoshi, and T. Watanabe, 
Nature 388, 431 (1997). 

8 R. Wang, K. Hashimoto, A. Fujishima, M. Chikuni, E. 
Kojima, A. Kitamura, M. Shimohigoshi, and T. Watanabe, 
Adv. Mat. 10, 135 (1998). 

9 K. Seki and M. Tachiya, J. Phys. Chem. B 108, 4806 
(2004). 

10 T. Thorsen, S.J. Maerkl, and S.R. Quake, Science 298, 580 
(2002). 

11 D. Juncker, H. Schmidt, U. Drechsler, H. Wolf, M. Wolf, 
B. Michel, N. de Rooij, and E. Delamarche, Anal. Chem. 
74, 6139 (2002). 

12 C. Bauer, S. Dietrich, and A.O. Parry, Europhys. Lett. 47, 
474 (1999). 

13 C. Bauer and S. Dietrich, Phys. Rev. E 60, 6919 (1999). 

14 M. Brinkmann and R. Lipowsky, J. Appl. Phys. 92, 4296 

(2002) . 

15 M. Schneemilch, N. Quirke, and J.R. Henderson, J. Chem. 
Phys. 118, 816 (2003). 

16 M. Sprenger, F. Schlesener, and S. Dietrich, Phys. Rev. E 
71, 056125 (2005). 

17 S. Kondrat, A. Poniewierski, and L. Harnau, Liquid Crys- 
tals 32, 95 (2005). 

18 Q. Wang, Macromol. Theory Simul. 14, 96 (2005). 

19 Y. Tsori and D. Andelman, Europhys. Lett. 53, 722 (2001). 

20 Y. Tsori and D. Andelman, Macromolecules 34, 27190 
(2001). 

21 Y. Tsori and D. Andelman, Interface Science 11, 259 

(2003) . 

22 D.J. Diestler, M. Schoen, J.E. Curry, and J.H. Cushman, 
J. Chem. Phys. 100, 9140 (1994). 

23 H. Bock and M. Schoen, J. Phys.: Condens. Matter 12, 
1545 (2000). 

24 H. Bock, D.J. Diestler, and M. Schoen, J. Phys.: Condens. 
Matter 13, 4697 (2001). 

25 H. Bock, D.J. Diestler, and M. Schoen, Phys. Rev. E 64, 
046124 (2001). 



2b H. Greberg and G.N. Patey, J. Chem. Phys. 114, 7182 

(2001) . 

27 S.D. Overduin and G.N. Patey, J. Chem. Phys. 117, 3391 

(2002) . 

28 H.B.G. Casimir, Proc. Koninkl. Ned. Akad. Wetenschap, 
BS1, 793 (1948). 

29 M.E. Fisher and P.-G. de Gennes, Compt. Rend. Acad. 
Sci. B 287, 207 (1978). 

30 M. Krech and S. Dietrich, Phys. Rev. Lett. 66, 345 (1991); 
67, 1055 (1991). 

31 M. Krech and S. Dietrich, Phys. Rev. A 46, 1886 (1992). 

32 M. Krech and S. Dietrich, Phys. Rev. A 46, 1922 (1992). 

33 M. Krech, Phys. Rev. E 56, 1642 (1997). 

34 F. Karimi Pour Haddadan, F. Schlesener, and S. Dietrich, 
Phys. Rev. E 70, 041701 (2004). 

35 HW. Diehl, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lebowitz (Academic, London, 
1986), Vol. 10, p. 75. 

36 K. Binder, in Phase Transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lebowitz (Academic, London, 
1983), Vol. 8, p. 1. 

37 T.C. Lubensky and M.H. Rubin, Phys. Rev. B 11, 4533 
(1975). 

38 Y. Okabe and K. Ohno, Phys. Rev. B 30, 6573 (1984). 

39 A.J. Bray and M.A. Moore, J. Phys. A 10, 1927 (1977). 

40 T.W. Burkhardt and H.W. Diehl, Phys. Rev. B 50, 3894 
(1994). 

41 J.O. Indekeu, M.P. Nightingale, and W.V. Wang, Phys. 
Rev. B 34, 330 (1986). 

42 E. Stanley, Phase Transitions and Critical Phenomena, 
Oxford University Press (1971). 

43 B.J. Schulz, K. Binder, and H. Miiller, Phys. Rev. E 71, 
046705 (2005). 

44 M.E. Fisher and H. Au-Yang, Physica A 101, 255 (1980) 

45 M.E. Fisher and H. Nakanishi, J. Chem. Phys. 75, 5857 
(1981) 

46 J.C. Collins, Phys. Rev. D 14, 1965 (1976). 

47 L.S. Brown, Ann. Phys. 126, 135 (1980). 

48 M. Fukuto, Y.F. Yano, and P.S. Pershan, Phys. Rev. Lett. 
94, 135702 (2005) 

49 T. Getta and S. Dietrich, Phys. Rev. E 47, 1856 (1993). 

50 R. Garcia and M. H. W. Chan, Phys. Rev. Lett. 83, 1187 
(1999). 

F. Karimi Pour Haddadan and S. Dietrich, e-print, 
cond-mat/0512217 the notations have to be translated ac- 
cording to: L -> d, S -> C/2, D -> 5, Ci -> X 2 , ( 2 Xi. 



